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Abstract. - The phase behavior of binary fluid mixtures of hard hyperspheres in four and 
five dimensions is investigated. Spinodal instability is found by using a recent and accurate 
prescription for the equation of state of the mixture that requires the equation of state of 
the single component fluid as input. The role played by the dimensionality on the possible 
metastability of the demixing transition with respect to a fluid-solid transition is discussed. The 
binodal curves in the pressure-chemical potential representation are seen to lie on a common 
line, independent of the size ratio. 



The phase diagram of the simplest model of a binary mixture, that of additive hard spheres, 
has not been completely determined up to now and remains as an interesting open problem. 
Ever since the solution of the Ornstein-Zernike equation with the Rogers- Young closure for 
such a mixture given by Biben and Hansen jjj, the until then widespread and unchallenged 
belief that these mixtures cannot phase separate into two fluid phases (which was supported by 
the analytical results of the Percus-Yevick approximation j|) has been seriously questioned. 
In fact, a wealth of papers dealing with this and related demixing problems have recently 
appeared in the literature but the issue is still far from being settled. One of the key 

aspects of this fluid- fluid transition in hard-sphere mixtures is that, if it exists at all, it must 
be entropy driven. In contrast, in other mixtures such as molecular mixtures, the transition 
is driven by an asymmetry between the like- and unlike-particle interactions, i.e. it is energy 
driven, and at least the qualitative phase behavior is rather well characterized Q . 

Entropy driven demixing was nicely shown by Frenkel and Louis to occur in a simple 
lattice model for a binary hard-core mixture ||. Experimental work has provided evidence 
that results for (asymmetric in size) hard-sphere mixtures approximate rather well the ther- 
modynamic properties of real colloidal systems, so it has been suggested that a likely physical 
mechanism behind the entropy driven demixing transition is osmotic depletion. This notion, 
first introduced to explain phase separation in polymer-colloid mixtures, is included in the 
scaled particle theory of Lekkerkerker and Stroobants [Q , which predicts a fluid-fluid spinodal 
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instability for highly asymmetric mixtures but whose domain of validity is not clear yet. The 
depletion picture is also present in other more recent studies |^,^]. 

A different approach to the demixing problem relies on the use of equations of state (EOS) 
for the hard-sphere mixture. It is well known that the popular Boubhk-Mansoori-Carnahan- 
Starling-Leland (BMCSL) EOS of state [|l0| does predict no demixing. On the other hand, 
Coussaert and Baus |J have proposed an alternative EOS with improved virial behavior that 
predicts a fluid-fluid transition but at such high pressures that it is argued to be metastable 
with respect to a stable fluid-solid transition. In this Letter we also take the EOS approach for 
binary mixtures of hard additive hyperspheres . By looking into the higher dimensionality, one 
can gain insight into any thermodynamic phenomenology that extends to such dimensionality 
p| . We will show that indeed these systems clearly exhibit a spinodal instability behavior. 

Let us consider an ./V-component mixture of (additive) hard spheres in d dimensions. The 
total number density is p, the set of molar fractions is {xi}, and the set of diameters is {ci}. 
The packing fraction is r/ = 53il=i ? 7i = v dP{ ljd )i where r\i — VdPidf is the partial packing 
fraction due to species i, pi — pxi is the partial number density corresponding to species 
i, Vd — (7r/4) d / 2 /L(l + d/2) is the volume of a d-dimensional sphere of unit diameter, and 
(<7 n ) = J2iLi x i (7 i' ■ I n previous work jl2| we have proposed a simple EOS for the mixture, 
Z m {fl)y consistent with a given EOS for a single component system, Z s (r]), where Z = p/ pksT 
is the compressibility factor, p being the pressure, T the absolute temperature and ks the 
Boltzmann constant. Our EOS reads 
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In the one-dimensional case, eq. (Q) yields the exact result Z m (r]) = Z s (tt). Further, for d = 2 
and 3, it proved to be very satisfactory when a reasonably accurate Z s (r]) was taken @,|l|. 
In all instances examined with different Z s (r)), however, no demixing is found for the resulting 
EOS for the mixture in these dimensionalities. 

One-component fluids of hard hyperspheres (d > 4) have attracted the attention of a 



number of researchers over the last twenty years |11,14-1S|]. Amongst other similarities, they 
share the property of hard disks and spheres of exhibiting a first order freezing transition. 
This transition occurs at a packing fraction r/f that, relative to the close-packing fraction 7y cp , 
decreases monotonically with increasing d. Perhaps the most accurate proposals to date for 
Z s (r]) in d = 4 and d = 5 are the semiempirical EOS proposed by Luban and Michels Jl5| : 



Z B {rj) = l + b 2 v 
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where b n are (reduced) virial coefficients defined by the series Z s (r}) = 1 + J2^=i bn+iV™- 
Equation (||) is consistent with the (known) exact first four virial coefficients, regardless of 
the choice of C(v)- Luban and Michels observed that the computer simulation data favor a 
linear approximation for £(77) and by a least-squares fit procedure found £(77) = 1.2973(59) — 
0.062(13)7j/?7 C p for d = 4 and CO?) = 1.074(16) + 0.163(45)77/^ for d = 5. Substitution of 
eq. (||), using the above linear fits for ((77), into eq. (|l|) produces the EOS of a mixture of 
hard additive hyperspheres in d = 4 and d = 5 dimensions, respectively. These EOS are 
the basis of our subsequent analysis. As an illustration of the accuracy of our EOS, in Fig. 
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Fig. 1 - Compressibility factor for three equimolar mixtures in 4D and 5D systems. Lines are the 
theoretical results and symbols are computer simulation data Jl9| |. 

[l] we display the comparison of the density dependence of Z m (r]) as obtained in this work 
with the results of computer simulation ]l9| ] for three equimolar binary mixtures and different 
values of the diameter ratio a = e^/ci, both in 4D and 5D. As clearly seen in the figure, the 
agreement is virtually perfect. This agreement has been amply confirmed for other values of 
the concentration x\ and the diameter ratio a ]l9| ]. 

We next consider the instability of a binary mixture (N = 2) by looking at the Helmholtz 
free energy per unit volume, a, which is given by 



where Xi is the thermal de Broglie wavelength of species i. We locate the spinodals through 
the condition a\\axi — a\i — 0, with = d 2 a/dpidpj. Due to the spinodal instability, the 
mixture separates into two phases of different composition. The coexistence conditions are 
determined through the equality of the pressure p and the two chemical potentials [i\ and 
H2 in both phases (/i.; = da/dpi), leading to binodal curves. It should be noted that in our 
calculations we have considered the range 771+772 < ?7cp (t7c P = 7r 2 /16 ~ 0.617 in 4D and 
7? C p = 7r 2 v / 2/30 ~ 0.465 in 5D |l5|) and that without loss of generality we have set fc^T = 1. 
Also, for future reference it is useful to recall that the freezing transition in the same systems 
occurs at a packing fraction rjf ~ 0.31 and rji ~ 0.20, respectively (T(f |. 

Figures || and || show the results for the spinodals and binodals in the 1J2 vs 771 plane and 
for the binodals in the paf vs 771 plane for different values of a and for d = 4 and d = 5, 
respectively. The location of the critical point tends to go down and to the right in the 1J2 
vs r\\ plane as a decreases. Also, in both dimensionalities and for a given a, the binodal 
and spinodal curves become practically indistinguishable to the right of the critical point. In 
contrast, while the binodals with the bigger a are always above those corresponding to a lower 
diameter ratio, to the left of the critical point the different spinodal curves cross each other 
for d = 4. The value of the critical pressure p c (in units of cr^ d ) is not a monotonic function 
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Fig. 2 - Spinodal curves (upper panel: lines) and binodal curves (upper panel: symbols; lower panel: 
lines) in a 4D system. The filled symbols are the critical consolute points. 

Fig. 3 - Same as Fig. El but for a 5D system. 



of a; its minimum value lies between a = 1/3 and a = 1/2 when d = 4, and it is around 
a = 3/5 for d = 5. This non-monotonic behavior was also observed for hard spheres [||. 

It is conceivable that the demixing transition in binary mixtures of hard hyperspheres 
in four and five dimensions described above may be metastable with respect to a fluid-solid 
transition (cf. the values of rji) , as it also occurred with the case of hard spheres || . In fact the 
value of the pressure at the freezing transition for the single component fluid is pia d ~ 12.7 
(d — 3), 11.5 (d = 4), and 12.2 (d — 5), i.e. pf<7 d does not change appreciably with the 
dimensionality but is clearly very small in comparison with the critical pressures p c <J d we 
obtain for the mixture; for instance, p c af ~ 1000 (d = 3, a = 3/10) 300 (d = 4, a = 1/3), 
and 123 (d = 5, a = 3/5). However, one should also bear in mind that, if the concentration 
X\ of the bigger spheres decreases, the value of the pressure at which the solid-fluid transition 
in the mixture occurs in 3D is also considerably increased with respect to pi (cf. Fig. 6 of 
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Ref. ||). Thus, for concentrations x± ~ 0.01 corresponding to the critical point of the fluid- 
fluid transition, the maximum pressure of the fluid phase greatly exceeds pf. If a similar 
trend with composition also holds in 4D and 5D, and given that the critical pressures become 
smaller as the dimensionality d is increased, it is not clear whether the competition between 
the fluid-solid and the fluid-fluid transitions in these dimensionalities will always be won by 
the former. The point clearly deserves further investigation. 

An interesting feature arising in our results, and which to our knowledge has not been 
discussed up to now, must be mentioned. There is a remarkable similarity between the binodal 
curves represented in the paf-iji and in the ^i~r\\ planes (not shown). By eliminating r\\ as 
if it were a parameter, one can represent the binodal curves in a /i^ vs paf plane. This is 
shown in Figs. [| and || where we have chosen the origin of the chemical potentials such as 
to make Xi = Oi. The binodals in the in—paf plane collapse into a single universal curve 
(which is in fact almost a straight line) for each dimensionality (d = 3, d = 4, and d = 5). 
We have obtained the data for d = 3 by using the EOS due to Coussaert and Baus ](|. For 
each dimensionality, the curve /ii = F\ (pcrf) is seen to be independent of the asymmetry in 
sizes as the contributions from each a to a portion of the curve overlap. Something similar 
occurs with the curve \ii = Fi (paf) and the calculations indicate that F\ = F2 in the 
sense that if plotted on the same plane, both curves seem to probe different regions of the 
plane but otherwise overlap reasonably well. This is illustrated by Fig. || for the case d = 4, 
where the points of ^,2 corresponding to a = 3/4 connect with those of n\ for a = 1/3 and 
a = 1/2. These properties could reveal the existence of a certain geometrical regularity that 
might characterize the entropy driven transitions. One could then state that all binary hard 
(hyper)sphere mixtures belong to the same universality class since, for fixed d, there are some 
properties such as the binodal curves in which the asymmetry is to a certain extent irrelevant. 
We also find an upper quasi-universality: for the range of common values, the universal curve 
in 3D lies close but somewhat above the universal curve in 4D, and in turn this latter lies 
above the one in 5D (note, however, that the origin of ^ is actually a matter of choice). All 
three curves approach straight lines of slopes 0.44 [d = 3), 0.45 (d = 4), and 0.43 (d = 5) in 
the region of their higher values. 

The results presented in this Letter confirm that the geometrical effects of osmotic depletion 
become more important as the dimensionality increases in much the same way that in 3D they 
are more important for parallel hard cubes than for hard spheres. It could be argued that 
the results we have presented so far (except for the universality, which is also found with a 
completely independent EOS in 3D) strongly depend on our prescription (U) for the EOS of the 
mixture and on the use of eq. (j^). To address the first point, we have carried out calculations 
using the van der Waals one-fluid theory with the Luban-Michels EOS in 4D and found no 
spinodal instability for r/i + 7/2 < i] cp . Concerning the second point, we have also performed 
calculations using Eq. (Q) for Z m (rj) with the EOS proposed by Baus and Colot [jl7| for the 
single component fluid. In this instance we also find spinodal curves which do not coincide but 
are qualitatively similar with the ones obtained using eq. (||) . The sensitivity of the spinodal 
curves with respect to the EOS was also noted in the case of hard spheres . Finally, we want 
to mention that our approach opens up the possibility of examining whether there is a critical 
dimension above which the fluid-fluid transition in hard hypersphere mixtures becomes stable, 
including the role that the asymmetry plays in it. Also we hope that the universal feature 
that we have described will stimulate other studies to clarify its physical origin and further 
implications. 
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Fig. 4 - Binodal curves in the fii vs pof plane for d — 3 (circles: 02/01 = 1/5, squares: oi/oi = 1/10), 
d — 4 (symbols as in Fig. |^), and d — 5 (symbols as in Fig. ^|). Note that for each size ratio the 
binodal curve is restricted to p > p c . 

Fig. 5 - Binodal curves in the /u, vs paf plane for a 4D system. The case i = 1 corresponds to 
Hi > 200 and is represented by the large symbols, while the case i — 2 corresponds to /ij < 200 and 
is represented by the small symbols. The inset is a magnification of the box enclosed by the dashed 
lines and shows [12 vs po\ for 02/01 = 1/2, 1/3 and 1/5. 
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